Exact out-of-equilibrium central spin dynamics from integrability 
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We consider a Gaudin magnet (central spin model) with generic time-dependent exchange cou- 
plings. We explicitly show that the Schrodinger equation is analytically solvable in terms of general- 
ized hypergeometric functions for particular choices of the time dependence of the coupling constants. 
Our method establishes a new link between this system and the 5*17(2) Wess-Zumino-Witten model, 
and sheds new light on the implications of integrability in out-of-equilibrium quantum physics. As 
an application, a driven four-spin system is studied in detail. 



Introduction - The problem of describing the coherent 
out-of-equilibrium evolution of driven many-body quan- 
tum systems has attracted a great deal of attention in 
recent years. This interest was spurred by the recent 
advances in cold atoms and semiconductor physics, 
which made experimental observations possible. The 
attention of the community has mostly been devoted 
to the investigation of two limiting cases: the quench 
regime, where the variation of the parameters of the 
Hamiltonian is very fast with respect to all the other 
time scales of the problem, and the adiabatic regime, 
where it is slow (see [l[ for a recent overview). Outside 
of these two extreme situations, very little is known, 
either analytically or numerically. This is unfortunate, 
because this driven regime has the most potential for 
novel physics. 

In this letter, we investigate a model of N spin— 1/2 
degrees of freedom coupled by time-dependent exchange 
parameters Ji(t), 



N-l 



Hcs(t) = J2 J ^) S ° ' S " 



(1) 



Label refers to the 'central spin' which is coupled to 
the N — 1 other (mutually uncoupled) spins. For time- 
independent couplings, this model is known as a Gaudin 
magnet/central spin model Q. Crucially, this Hamilto- 
nian is directly relevant to experiments in quantum dots 
[H, 0] and Nitrogen vacancy (NV) centers in diamond [j| , 
in which time-dependent couplings are intrinsic to the 
experimental protocols (respectively via time-dependent 
gate voltages or external electromagnetic fields). 

The aim of this Letter is twofold. On the one hand, we 
identify a time-dependent protocol for which it is possi- 
ble to obtain analytical information (i.e. the exact many- 
body wavefunction) for (TTJ). While the requirements of 
this protocol are restrictive, they nonetheless allow to go 
beyond the adiabatic or sudden approximation. On the 
other hand, our technique points to an intriguing link be- 
tween the time dependent central spin Hamiltonian and 
the Wess-Zumino-Witten (WZW) model, a well-known 



conformal field theory (CFT), opening the door to fur- 
ther applications of CFT techniques to driven nonequi- 
librium physics. 

Main Results - Our starting point is the fact that the 
conformal blocks of the Wess-Zumino-Witten model @ 
satisfy the Knizhnik-Zamolodchikov equations Q- For a 
SU(2)k WZW model, these equations read 



Zi Zj 



^n{zq, ■ . -,Z N -l) 



= 0, 



(2) 



where ^ N (z , . . . ,z N -i) = (<p(z ) . . . ip(z N -i)) is the 
N-point holomorphic conformal block of primary field tp, 
while A; is a number known as the level of the Kac-Moody 
algebra. If k is a positive integer, the WZW model 
is a rational CFT. Interestingly, there exist integral 
representations of solutions to the KZ equations that can 
be analytically continued to any nonzero complex k [8|,|9|. 

Choosing k = iv where v € 1R and considering the 
ansatz ipN(t) — ^ N(zo(t), ■ ■ ■ > z N-i{t)) for a many- 
body wavefunction, we see that ipN(t) can in fact be 
reinterpreted as a time-dependent Schrodinger equation 
itv4>M{t) — H(t)ipN(t) with Hamiltonian 



H(t) 



EE 

» 3^ 



hzjjt) Sj-Sj 

v Zi(t)-Zj(t)' 



(3) 



Therefore, if zq (t) is chosen to be the sole time-dependent 
parameter, ipN(t) — ^N(zo(t), Z\, . . . , -zat-i) solves the 
time-dependent Schrodinger equation for Hamiltonian 
([1]) with couplings 



Ji(t) 



hz (t) 



v {z (t) - Zi) ' 



(4) 



It is shown in the Supplementary Information (SI) that 
this choice of time-dependent parameters Zi(t) is uniquely 
dictated by the form of ([1]). Notice that the hermiticity 
of the Hamiltonian forces all the Zi to be on the same line 
in the complex plane (for example, we can take them to 
be all real). 
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Let us emphasize the main features of our approach. 
First of all, the main ingredient for an explicit solution of 
the time-dependent Schrodinger equation is the integral 
representation for the solutions of the KZ equations. 
These representations rely crucially on the integrability 
of the time-independent Hamiltonian, i. e. the off shell 
Bethe equations. Therefore, the solubility of the time- 
dependent Schrodinger equation seems to be a signature 
of the underlying integrability of the model that survives 
also when the coupling are time dependent. Indeed, this 
interpretation is confirmed by the fact that, as we will 
discuss later on, the solubility of the time-dependent 
Schrodinger equation is not a special feature of the cen- 
tral spin model: our arguments apply also to the broader 
class of XXZ Gaudin magnets. Moreover, our results 
establish a new connection between the SU(2) WZW 
model (admittedly, for the quite unusual imaginary k 
case) and the time-dependent central spin Hamilto- 
nian. It is worth to note here that the SU(2) WZW 
model is known to be related to integrable[10f and non 
integrable[ll|, |l2j time-independent spin Hamiltonians. 
Finally, it is important to stress that-by construction- 
this approach works only if the time dependence of the 
Jiit) is finely tuned: essentially, the time evolution is 
"geometric" , i.e. J dtH(t) can be written as a curvilinear 



integral in the space of the Zj's. 

A simple example - Before presenting the details of our 
general solution, we would like to focus our attention on 
the very simple case of a central spin Hamiltonian with 
four constituents, 



H(t) 



1 3 

i} £ — ' 



i *>(*) 



-So • Sj 



(5) 



Since the total spin is conserved by the time evolution, 
we can restrict ourselves to the subspace with constant 
S 2 = (J^i Si) 2 . Here, we would like to show that, indeed, 
the WZW correlators provide solutions that describe the 
whole zero spin subspace. 

The computation of the four point conformal blocks 
^(zo, zi,z 2 ,z 3 ) of the WZW model is a standard exer- 
cise of CFT (see Q). The detailed calculation is reported 
in the SI, where ^^(zo, Z\,z%, 23) is expressed in terms of 
the standard hypergeometric functions 2F1 (a, b, c, x) [HI . 
We introduce the parametrization ^(^Oj Z\, Z2, Z3) = 
[(z — z 3 )(zi — Z2)]" f (x) , where /(x) is a function 
of the anharmonic ratio x = ( z °~ Zl )( Z2 ~ Z3 ) anc [ expand 

f(x) = [x(l — x)]~3T Y^i=i Gi{x)\vi) on a basis of the 
S 2 = subspace given by the two states 



V2 



V2 



V2 



(6) 



1 

V3 



+ -) + !-+) 

V2 



V2 



We thus can write Gi(x) = Yli=i <HVJi(x), where ci,2 are 
constants determined by the initial conditions, while 

^) = ^(-|>-^>nH' (7) 

/ 1 2k 1 \ 
Wi(x) = (~x)— 2 Fi f 1 - — , 1 + — + -,xj , 

and G 2 (x) = [3Gi(x) +4kxG[(x)}. There- 

fore, as discussed above, the wavefunction ^64 (i) = 
^4(zo(t), Zi,Z2, 23) (with k — iv) is a solution of the time- 
dependent Schrodinger equation for Hamiltonian ([S]). 

As an example, let us consider the following proto- 
col. At time t=0, the spins Sj, j = 1,2,3, are at a dis- 
tance j from the central spin So- Their couplings Jj are 
taken to be proportional to j~ 3 (dipolar interaction) or 
to exp (— j 2 ) (shell model). Subsequently for t > 0, the 
coupling constants decrease inverse linearly in time. We 
can thus model this situation with zq = tut and Zj = —j 3 
(dipolar interaction) or Zj = — exp (j 2 ) (shell model), 
.7 = 1,2,3. 



The first thing to analyze is the completeness of the 
solution, i.e. if the space spanned by the conformal 
blocks solution is bidimensional. It is shown in the 
SI that the absolute value of the determinant of the 
matrix My(x(t)) such that Gi(x(t)) — J2j Mij{x(t))cj 
is actually constant and nonvanishing for t G [0, +00), 
thus proving that this family of solutions spans the 
whole subspace of zero total spin. This fact is unrelated 
to our choice of zo (t) and Zi and remains true for 
any parametrization (see SI). As an application, an 
interesting quantity to look at is the modulus square of 
the overlaps of the wavefunction with the basis vectors 
\vi) , i.e. di(t) = \ {vi\ip(t))\ 2 which are simply computed. 
Two interesting examples are shown in Fig. ([T]) for the 
dipolar interaction (left) and for the shell model (right). 
In both cases, the initial condition is chosen in such a 
way that it is possible to observe one (dipole interaction) 
or two (shell model) crossings of the overlaps a, (t) . 

General Solution- It this section, we would like to sum- 
marize our strategy for getting an integral representation 
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FIG. 1. Time evolution of the modulus square of the overlaps ai(£) (solid line) and a,2(t) (dashed line) of the wavefunction 
with the basis vectors for the dipole interaction (left) and for the shell model (right). In these plots uj — 10. In the inset: a 
magnification of the shell model plot for t £ [0, 10]. 



of the solution to the KZ equations ^ or- more precisely- 
a generalized version of these equations. Our arguments 

I 



are very close to the ones in [8|, Il4|, [l_5j . The XXZ Gaudin 



magnets are defined from the Gaudin algebra 



[S»(Ai), S Z (X 2 )} - i [Z(X 1 - X 2 )S X (X 1 ) - X(Ar - X 2 )S X (X 2 )} , 
[5 z (Ar), S X (X 2 )} = i [X(Ai - X 2 )S y {Xi) - Z{X X - X 2 )S y {X 2 )\ , 



[s x (x 1 ),sy(x 2 )]=ix(x 1 

[S a (X 1 ),S a (X 2 )} = 0. 



A 2 )[5 Z (A! 



S Z (X 2 )}, 
(8) 



Here, X and Z are odd functions, with Res [X(z)] z _ — 
Res [Z(z)] z=0 = g, while the A^'s are complex numbers. 
Notice that X and Z are not arbitrarily functions, but 
they have to satisfy a set of quadratic equations that 
come from the Jacobi identities for the generators of 
the Gaudin algebra (jHJl. The solutions to these equa- 
tions are known, and the simplest one is the rational one 
X(Xi,X 2 ) = Z(Xi — X 2 ) = Xi 9 _ X2 (for a detailed discus- 
sion of Gaudin magnets, the reader is referred to [16|). 
For example, in order to describe a spin or fermionic sys- 
tem, the su(2) representation is useful (S ± — S x ± iS v ) 



S ± (X,z) 



]X(\-Zi)S* 



S z (X,z) 



Z{X- Zi )S*. 



(9) 



(10) 



Here Zi are a set of complex numbers {the disorder vari- 
ables) that are directly linked to the coupling constants 
of the Hamiltonian, while Si are the familiar spin opera- 
tors. Instead, a bosonic system is described by a su(l, 1) 
representation (z = Zo, . . . Zn-%) 

S ± (X,z) = ±Y,X{X-z t )K±, (11) 



S*(X,z) 



1 



)Kf 



(12) 



where the KiS satisfy a su(l, 1) algebra. Of course, 
mixed representations are also available, in order to de- 



r 



scribe a system where the spin degrees of freedom interact 
with a bosonic bath (i.e. the Dicke model and its gen- 
eralizations). In the following, we will use the notation 
Af L = Sf /Kf respectively for the su(2) representation or 
the su(l, 1) one. 

From the Gaudin algebra, it is possible to define the 
generating function 



S a (uj,z)S a (uj 1 z), 



and the integrals of motion 



(13) 



(14) 



These H^s are a family of commuting operators, and 
each of them can be considered as the Hamiltonian of 
a quantum system. Since H(w, z) can be diagonalizcd 
using the algebraic Bethe Ansatz, the Hamitonians Hi 
are exactly solvable. In the rational case, the integrals of 
motion for a spin system reduce to 



2ff 



(15) 



that is a central spin Hamiltonian with a Zeeman term 
for the central spin: when the magnetic field is zero it 
reduces to the operators appearing in the KZ equations 
@. The main statement of this section is that, for the 



4 



spin as well as for the bosonic representation, it is pos- 
sible to construct an explicit solution of the generalized 
KZ equations 



k^{z)=H i {z)* N {z). 



(16) 



This solution, as we have explained in the first section 
of this Letter, is not simply a mathematical curiosity, 
but it can indeed be used to solve a time dependent 
Schrodinger equation: ipN(t) = ^N(zo(t), ■ ■ ■ , zjsr-i(t)) 
satisfy the Schrodinger equation with Hamiltonian 
H(t) = jEj^W^'tiW)- Quite nicely, this solution is 
based on the integrability of the model. Let us introduce 
the Bethe state (A = Ai, . . . , Am ) 



M 



l*(A^)> = l[S+(\ a ,z)\0), 



(17) 



where the vacuum |0) is annihilated by S~(X) and is an 
eigenstate of S Z (A) (Af |0) = di\0)). It can be proven 
that this state obeys the off-shell Bethe equation 



Hi(z)\$(\,z)) = h t (X,z)\<£(X,z)) + (18) 
+ X>(A«-^)/«(A,z)A+ [] S+(Xp,z)\0). 



The functions hi and f a can be derived from a Yang- 
Yang action: more precisely, hi = 
S is defined as 



fa = where 



M 



EE^ 



Aq) 



M A 



M M 

EE 

a /3^a 



T{\/3 - Aq) 

2.9 : 



(19) 



where T(u) = J™ dzZ(z). Usually, one imposes the on- 
shell condition /" — (Bethe equations), thus obtaining 
a basis of eigenstates of the Hamiltonian. Here instead, 
we take advantage of the existence of the action S and 
we define 



dX e £ % ii |$(A,z)), 



(20) 



where the closed contour 7 is chosen in such a way 
that the branch of the integrand at the end point 
of 7 is the same as that at the initial point. It is 
quite easy to show that (|2"0")) is indeed a solution of 
(IT6| . Notice that due to the multi-valuedness of the 
integrand, the path of integration is usually highly 
nontrivial, this being the major technical difficulty 
of our approach. In the rational case, these integrals 



represent multivariable hypergeometric functions 13, 



The k — > limit- From ([2H)l . it is quite clear that in 
the k — > limit, the only contribution to the integral 
comes from the stationary points of the action (TT9")) , i.e. 
from the on-shcll Bethe state Therefore, it is quite 
interesting to discuss the physical meaning of this limit 
for our time-dependent Schrodinger equation. The most 
natural interpretation of this limit is as an adiabatic one. 
As an example, let us consider the central spin limit with 
coupling constant (jj). If we parametrize k = iv and 
Zo(t) = F(p,Qvt), where F is an arbitrary function, we 



have a central spin model with coupling constants 



j l (t) = hn 



F'(n vt) 
F(n vt)- Zt 



(21) 



The time scale of Ji(t) is ft^v, and so when v — » 
the system is in an adiabatic regime. Notice that, in- 
deed, the contribution from the stationary points of (|2U)) 
agrees completely with the usual quantum adiabatic the- 
orem: the stationary condition f a — J^- = imposes the 
Bethe equations, thus selecting the instantaneous eigen- 



state of the Hamiltonian, while exp 



k 



/a=0 



is the 



corresponding dynamical phase. 
Quite interestingly, the k — > limit can be inter- 
preted also as a semiclassical limit, if we use a different 
parametrization. Indeed, if k = i-j^i where Aq has the 
dimensions of an action, we have a central spin model 
with coupling constants 



Ji{t) = A - 



zo{t) 



zo(t) - Zi ' 
and k — > is equivalent to Ao h. 



(22) 



Conclusions- In this Letter, we have studied a class 
of time-dependent Hamiltonians which possess many- 
body wavefunctions given by solutions to the associ- 
ated Knizhnik-Zamolodchikov equations. The underly- 
ing time-independent integrability and link with CFT al- 
low us to provide an explicit integral representation of the 
solution to the time-dependent Schrodinger equation. It 
is worth to notice that, while these integrals representing 



multivariable hypergeometric functions are quite compli- 
cated, they have been extensively studied in the mathe- 
matical literature. For a small number of particles, these 
reduce to the familiar hypergeometric functions, allow- 
ing us to easily study the dynamics of the system. While 
our construction works only if the time dependence of 
the Hamiltonian is finely tuned, it provides an intrigu- 
ing starting point for understanding the consequences of 
quantum integrability in time-dependent physics. 
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Supplementary Information 



THE CENTRAL SPIN MODEL AND THE WZW CFT 

In the text, we have argued that the conformal block of the SU (2)& ^n{zo, • • • , zn-i) can be used to construct the 
wave function tpff{t) = ^N{zo(t), Zi, ■ ■ ■ , Zn-i), that is a solution of the Schrodinger equation for the time dependent 
central spin Hamiltonian (k = iv) 



N-l 



Hcs{t) = £ 



hz (t) 



-So • Si. 



-f x v (z (t) - Zi(t)) 

This result follows easily from the K-Z equations satisfied by the conformal blocks 



(23) 



By.i 7.1 — 7.1 



$> n (z ,...,z N -i) =0. 



(24) 



Here, we would like to elaborate more on this point. In particular, one could wonder if there exists a more general wave- 
function ipN(t) — ^N(zo(t),zi(t), . . . , ZN-i(t)), where all the Zi{t) are time dependent, that satisfies the Schrodinger 
equation for a central spin model. The answer is essentially no. Indeed, ipN(t) — ^N(zo(t),zi(t), . . . , ZN-i(t)), satisfies 
a Schrodinger equation with Hamiltonian 



ffW = EEc, 3 (t)s r s 31 



Cij(t) 



v[Zi(t) - Zj{t)]' 



(25) 
(26) 



So, if we want to cancel out the couplings between the spins of the bath, we need to have z\{t) — z 2 (£)... = 
£jv-i(i)-This is equivalent to have the point zo(t) moving in time while zi, . . . , zn-i are fixed, once we take into 
account the invariance of the conformal block ^n(zo,zi,. . . ,zn-i) under global translations. 



THE FOUR SPINS CONFORMAL BLOCK 



In this section, we would like to derive explicitly the four point conformal block ^(^o, Z\, Z2, 23). The S 2 = 
subspace is spanned by two vectors 



+ -}-|-+>WI + -}-|-+> 



V2 



V2 



1 

7! 



+ --) + !-- + +) 



l + -> + l-+>WI + -> + l-+> 



V2 



V2 



It turns out that 



So - Si|«i> = So ■ Si|«2> = ||«2> 



So • S 2 |ui) = -^K) S -S 2 |u 2 ) 
Moreover, in this subspace S a \ip) — Sj\ip) — 0, hence 

- + S • Si + S • S 2 + S • S 3 = (S 2 = 0). 



(27) 



(28) 
(29) 

(30) 



So, let us consider the zq K-Z equation 

d So • Si So • S 2 Sq • S3 



8zq Zq — Z\ Zq — Z2 Zq — Z3 



^ i (z ,z ll z 2l z 3 ) = 0. 



(31) 
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The sum rule (|30[) teaches us that we can eliminate the dependence of the Hamiltonian on So • S3. Moreover, if we 
make the substitution 



* 4 (zo, Zi, Z 2 , Z 3 ) = [(^0 - Z 3 )(zi - ^2)] 4 M-ZO, Zl, «2, 23)}, 



the K-Z equation become 



.9 / 1 1 
1- h S • Si 

dz V zo - z 3 zo - 



+ Sq • S2 



zo - Z3 z - z 2 



\<p) = 0. 



Let us now introduce the new variable 



so 



Moreover, the following identity holds: 



(zo - zi)(z 2 - z 3 ) 
(zo - z 3 )(z 2 - zi)' 



<9z n 



zo — Zi z — z 3 



d_ 

dx' 



l-x 



z - zi zo - z 3 



zo - z 3 z - z 2 



(32) 



(33) 



(34) 



(35) 



(36) 



This identity can be proved by a brute force calculation (it's just a trivial algebraic calculation), but it is possible to 
obtain it in a clever way. First of all, we notice that the l.h.s. has poles in zo = z 3 and zo = z 2 (since x = 1), so we 
can deduce that 



l.h.s. 



Q(z ,zi,z 2 ,z 3 ) 
(z - z 2 )(z - z 3 ) ' 



(37) 



where Q is a polynomial. However, the l.h.s is zero only for z 2 = z 3 , hence Q = A(z 3 — z 2 ). Taking the limit z —¥ z 3 , 
we get A=l. 

Therefore, the K-Z equation reduces to 



d So • Si So • S 2 

dx x x — 1 



\ip(x,z 1 ,z 2 ,z 3 )) = 0. 



(38) 



Now, we can expand \tp) on our basis 



\ip(x,z 1 ,z 2 ,z 3 )) = y^Fj(x, zi,Z2,z 3 )\vi), 



»=i 



obtaining a system of differential equations for F± and i 7 ^ 



Ax 4 x - 1 
lFs , 1 



4 1 x — 1 



Fl -|- — Fn 

4 2 



(39) 

(40) 
(41) 



where the prime denote a derivative respect to x. Our aim now is to show that these equations reduce to a hyperge- 
ometric one. From the first one we get 



F a = -(*-!) 



yfiFi Ai 
1 — 1=^1 



and, the second equation reduces, after some algebra, to 

a(aj)J^'(a:) + (3(x)F[{x) + »/(x)Fi(x) = 0, 



(42) 



(43) 
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where 



a(x) 
0(x) 



21 + 4<i-) 



Xy/3 v3 
7(x) = v3 





1" 


1 1 


{ X 2 




4 x 



4 a; Ax -I 



or, equivalently, 



x(l - x)F['(x) 



+ (1 - i)« 



16 



x- 1 



F 1 (x)=0. 



This equation is indeed quite similar to the hypergeometric equation 

x(l - x)w"{x) + [c - (a + b + l)x] w'{x) - ab w(x) = 0, 



(44) 
(45) 
(46) 

(47) 
(48) 



but not identical, since the coefficient of the last term depends on x.How can we get rid of those terms?Let us define 
a function v{x) such that 



Fi(x) = x r v(x), 

where r is a complex number. Clearly, 

F[{x) — rx r ~ 1 v(x) + x r v'(x), 

Fi'{x) = r(r - l)x r ~ 2 v(x) + 2rx r - 1 v'(x) + x r v"(x), 
so, if we substitute these expressions in (I4T1) we get that the term proportional to v(x) is 



i 3(4i + l) 

r(r - 1) + -r 

V ' 2 16 



(49) 



(50) 
(51) 



(52) 



where we omitted terms regular for x = O.So, if we choose r = ^orr = l- |, we can eliminate the nasty - term. 
A similar argument applies for x — 1. So, let us define a function g(x) such that 

Fi{x) =x^(x - l)-ig(x). (53) 
It turns out that g(x) satisfy a hypergeometric equation 



x{l - x)g"(x) + (* - x)g'{x) - -g(x) = 0, 



(54) 



with parameters a 
equation, we have 



-b = |, c = i.So, if wi and u>2 are two linearly independent solutions of the hypergeometric 



Fi (x) = x * (x — 1) 4 [ciwi(x) + £2^2(2;)] 
where c\ and C2 are arbitrarily constant, determined by the initial condition. 



(55) 



A SIMPLE EXAMPLE 



Let us discuss now a simple example. At time t=0, the spins S.;, i — 1,2,3, are at a distance j from the central 



spin So. Therefore, their coupling with the central spin J{ is proportional to exp (- 
constant decreases linearly in time. So, we can model this situation if zq — ujt, and Zi 



. Then, for t > 0, the coupling 
exp (i 2 ) i = 1,2,3. Therefore, 



H(t) = LU 



So • Si Sq • S2 So • S3 



ojt + e cut + e 4 



(56) 
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FIG. 2. |det[M(a;)]| 



while 



x(t) 



r e 9 



Lot - 



Lot - 



-00,0). 



(57) 



So, there are 24 solutions of the hypergeometric equation in the complex plane [18j j . These solutions are characterized 
by different analytical properties: for example the hypergeometric function F(a, &, c, z) is analytic in 0, while F(a, b, a+ 
b + 1 — c, 1 — x) is analytic in 1. Of course, if different solutions are well defined in the same region, at most two 
of them can be linearly independent, since the hypergeometric equation is of second order. For our aims, a good 
choice of two linearly independent solutions is F(a, b, c, z) and z 1 ~ c F(l + a — c, 1 + b — c, 2 — c, z), since they have no 
singularity on the negative real line. So, we put 



w 1 (x) 

U>2 (x) = X 

while F2{x) = c±yi(x) + 022/2(2) is thus given by 

F 




+(.T-l)F(l--,l + -,l + l ,X 



•!.H r\ I -1,1- ^,2 



(58) 
(59) 



i,x ) + 



(60) 



and so 



Fi(x) = M ij (x)c j , M(x) 



wi(x) 

yi{x) 



w 2 (x) 



(61) 



As we have discussed in the text, the physics of this simple model could be quite interesting, with a double crossing of 
the overlaps \Fi\ 2 for some initial conditions. Moreover, as it is shown in fig. [2], the absolute value of the determinant 
of M(x) is constant for x < and non vanishing. This indeed imply that this class of solution spans the whole S 2 = 
subspace. More precisely, | det[M(x)]| is indeed a (non vanishing) constant in the connected domains (— 00, 0) , (0, 1) 
and (1, +00), with jumps at the singular points of wi and W2- The fact that, inside these domains, | det[M(x)]| is a 
non zero constant reflects the unitarity of the quantum time evolution, that implies 



|Fi(x)| 2 + |F2(x)| 2 = non vanishing constant, 



(62) 



or, equivalently 



M\x)M(x) = non vanishing constant, 



(63) 
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and therefore 

| det(M(x)| 2 = non vanishing constant. 



